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We consider tlie asymmetric simple exclusion process (ASEP) on a semi-infinite chain which is 
coupled at the end to a reservoir with a particle density that changes periodically in time. It is 
shown that the density profile assumes a time-periodic sawtooth-like shape. This shape does not 
depend on initial conditions and is found analytically in the hydrodynamic limit. In a finite system, 
the stationary state is shown to be governed by effective boundary densities and the extremal fiux 
principle. Effective boundary densities are determined numerically via Monte Carlo simulations and 
compared with those given by mean field approach and numerical integration of the hydrodynamic 
limit equation which is the Burgers equation. Our results extend straightforwardly beyond the 
ASEP to a wide class of driven diffusive systems with one conserved particle species. 
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I. INTRODUCTION 

Systems of driven diffusing particles attract attention 
because, despite their relative simplicity, they embrace 
a whole range of critical phenomena far from thermal 
equilibrium |l], [3, [!]■ One of the remarkable features 
of these systems is the appearance of phase transitions 
induced by spatial boundaries of an open system which 
exchanges particles with external reservoirs 0, @ . A clas- 
sical model where this can be studied in great detail is the 
so-called Asymmetric Simple Exclusion Process (ASEP) 
with open boundaries. This model describes the single- 
file random motion of particles with hard core exclusion 
and drift inside a finite system at the ends of which par- 
ticles can be extracted or injected with some rates. This 
model was first introduced for describing the kinetics of 
protein synthesis [1, 0| and has since then been general- 
ized in many ways for describing the motion of various 
kinds of molecular motors [1, H, . Due to its concep- 
tual simplicity it also plays a fundamental role in traffic 
flow theory [ill . [T^ and many other settings where driven 
diffusion of interacting particles plays a role. 

By now, the dynamics of the ASEP (as well as its sta- 
tionary bulk behavior) is rather well understood. For 
our purposes we note that the exact stationary distri- 
bution has been determined analytically [l^, [Tj], while 
the coarse-grained dynamics of shocks and localized ex- 
citations in the evolution of the particle density can be 
understood using hydrodynamic limit equations [l^ . 
The latter provides a full description of the evolution 
of the local density under Eulerian scaling. It has been 
shown rigorously [17,] to be given by the famous Burg- 
ers equation used for the description of the dynamics of 
shocks in dissipative systems [185 . 

The vast body of knowledge about the ASEP has 
been obtained for time-homogeneous conditions where 
the boundary rates are kept constant in time. In con- 
trast, very little is known when the environment of this 
open system changes non-adiabatically in time on scales 
that are comparable to the macroscopic Eulerian hy- 



drodynamic regime. This has to be modelled by time- 
dependent boundary rates which, to our knowledge, has 
not yet been attempted for the ASEP with open bound- 
aries. It is the purpose of this work to report simulation 
results for a natural time-periodic setting and to analyze 
these data in the framework of the hydrodynamic theory. 

The paper is organized as follows. In Sec. |TT]we define 
the model and present our simulation data. The data for 
the Eulerian low-frequency regime of a semi-infinite sys- 
tem are then analyzed and explained by extending hydro- 
dynamic theory to incorporate time-dependent boundary 
conditions (|TTT}. This analysis allows us to predict the 
phase diagram of an open finite system with two bound- 
aries. This prediction and its numerical verification is 
given in Sec. IIVI We end with a summary of our results 
and some conclusions dVl). 



II. ASEP MODEL ON A SEMILINE WITH 
TIME PERIODIC BOUNDARY 

We consider the ASEP defined on a semi-infinite chain 
{fc g Z, , /c < 0} with a right boundary site fc = 0. Each 
site of the chain is either occupied by one particle or 
empty. We denote the local occupation number by rik € 
{0, 1}. Particles attempt to jump to the right or to the 
left neighboring site after an exponentially distributed 
random time with parameter p -I- g, normalized asp + q = 
1. The rate at which a particle attempts to hop to the 
right (left) isp (q). If the target site is empty, the attempt 
is successful and the particles moves. Otherwise, it does 
not jump (hardcore exclusion rule). At the boundary 
site fe = a particle can be extracted with the rate /3 
(if the boundary site is occupied) or be injected into it 
(if the boundary site is empty) with rate S. We choose 
f3 = p{l — pji) , S — qpB. so that the boundary may be 
thought of as being coupled to a reservoir of density pR 
p^ . We consider the boundary reservoir density pu to 
be a periodic function of time with frequency u, switching 
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Figure 1: Snapshots of density profiles averaged over 3*10^ 
histories, during red light t = 0.75T, t = r(lines) and during 
the green light t = I.IT, 1.2r, 1.3T, 1.5T (points). Parame- 
ters: T — 1000. Initial state is a homogeneous one with the 
density (nfc(O)) = 0.15. 



between the values pn, = and pji — I according to 

PR{t)^^{l + sign[sm{2Tn^t)]). (1) 

In traffic flow problems such a boundary condition mod- 
els a traffic light with pR{t) = during the "green light" 
half-periods t C [0, [T, and PR{t) = 1 during 

"red light" half-periods r^ed C [^,T], [^,2T], ... , where 
T = 1/v. In analogy with this we refer to Eq.([T]) as traffic 
light boundary condition. For a study of such a switching 
in a related system, see . 

To investigate the effects of this boundary condition 
on the dynamics of the ASEP on a semifine we have 
carried out Monte-Carlo simulations of the model for 
various frequencies v. We have concentrated mainly on 
the case of a totally asymmetric simple exclusion pro- 
cess (TASEP), g = 0,]3 = 1, and focused our interest 
on stationary behavior in the sense, that all macroscopic 
quantities behave periodically in time with period T . The 
initial state was prepared in an ensemble of particles ran- 
domly distributed with the density A, which is station- 
ary in the infinite system. In presence of periodically 
varying boundary rates, rather complicated dynamics is 
observed. For small ingoing fiuxes jin = A(l — A), the 
system develops jam at the boundary during the red light 
periods, which is dissolving completely during the green 
light periods, FigU] However, if A (and consequently, 
the infiow flux) exceeds some critical value A > Xc{v), 
the jam at the boundary is not dissolved completely, but 
starts to propagate inside the system. The amplitude of 
the shock front is not constant, it increases and drops 
during each red-green period, see FiglS] Consequently, 
shock front is not propagating steadily, but its velocity 
changes, and in particular it may advance and retract 
during each red-green period. The net shock advance 
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Figure 2: Snapshots of density profiles at times t' = t ~ = 
0, 50, 75, 100, at the coexistence line between LD and HD 
states jin — jout- Parameters: T = 100. The initial con- 
figuration of the system is a domain wall 0.2,0.7 positioned 
at k/T — —2 ( dotted piecewise straight line). The system 
was equilibrated for to = 5T — 500 before the measurements, 
and the averaging over 10^ histories was done. Thick line 
shows the density profile after t = 800 in the system prepared 
initially in HD state (taken from Fig[3}. Thin line shows the 
density profile averaged over several periods and many histo- 
ries after initial transient period of 5T and it looks (apart from 
boundary layer) like unbiased domain wall (0.2, 0.8 = p'jl'^ (u) 
in TASEP with constant boundary rates 



after a period T is determined by the mass conserva- 
tion, i.e. difference between the ingoing jinT and out- 
going fluxes of particles joutT, see Figs l3l4l While jin is 
a control parameter jin = A(l — A), jout is not and it 
is measured as the time-averaged particle flux through 
the boundary, jout depends on the frequency of traffic 
light ly and it is associated with the effective right bound- 
ary density p'r ^ {v) through the TASEP current-density 

relation jout ~ p'r/ {^) ^ p'r^ {^)^ as discussed be- 
low. If ingoing and outgoing fluxes are equal, the shock 
only "breathes" around its initial position see Fig[2] The 
density profile behind the shock front develops approxi- 
mately equidistant sawteeth-like structures with the de- 
caying amplitudes, see FigO The sawteeth profile is 
changing with time, but it regains its shape after each 
complete period, so that the shape depends only on the 
phase 7 C [0,r]. Apart from this phase dependence, 
illustrated e.g. on FigUl the sawteeth structure behind 
the shock front depends on frequency v (the rescaled saw- 
teeth become sharper with decreasing v ) but not on A, 
see Figs l2l3lil The particle density at some distance from 
the boundary approaches a constant value, p'^^^ {v). This 
value can be identified with an effective right bound- 
ary density in the following sense: the particle system 
(ASEP) behaves as if it was joined at the origin with 
the reservoir of particles with the density pR = iy) . 
As a function of frequency the averaged stationary flux 
through the boundary changes monotonically, which im- 
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Figure 3: Snapshots of averaged density profiles p{x,t) at 
t = Q (homogeneous state with the density (rife) = 0.5) and 
then at equal intervals of time (equal to two periods T) ti = 
200,400,600,800. Parameters: v = T"^ = 0.01, system size 
400, left boundary density 0.5, averaging is done over 3*10^ 
histories. Inset shows the snapshots of density profiles (rife) 
at times t = ti in case when instead of traffic light conditions, 
time-independent boundary conditions are applied with the 
same effective right boundary density pn = pj/^(i^) ~ 0.8. 



Figure 5: The averaged stationary density (for parameters 
chosen equal to effective right boundary density pj/^) in the 
TASEP versus the frequency of "traffic light" switch at the 
right boundary, v = , from Monte Carlo simulations. The 
parameters are: system size 400, the left boundary density is 
higher than 1/2, and the averaging is done over 30 histories 
and over 2 * 10^ Monte Carlo steps, after the equilibration. 
Points show the results of density averages, and the broken 
fine reports estimates from stationary flux measurements. 
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Figure 4: Same as in Fig[3] but for an initial homogeneous 
state with density p{x, 0) = 0.4. Dotted lines show the result 
on numerical integration of mean field equations (|6} with p — 
1. 



plies a monotonic change of Pj^ {v) ■ In Fig[5]we report 
the numerical results for p'jj/{i^) from the Monte-Carlo 
simulations. 



III. HYDRODYNAMIC THEORY AND 
SAWTOOTH STRUCTURE 

A. Limiting cases 

In order to understand the origin of these observations 
we first argue that if A > 1/2 then the right bound- 
ary reservoir density p^'^ (i^) must indeed be equal or 
larger than 1/2 for any value of z/, as observed: Dur- 
ing the "green light" periods /3 — 1 at most 1 parti- 
cle per time unit can exit, while during the "red light" 
periods /3 = no particles can exit. For frequencies 
large compared to the mean attempt rate for particle 
jumps (which is 1), a particle at the boundary site "sees" 
the reservoir with equal probability empty or fully oc- 
cupied, irrespectively of how long it has already stayed 
at the boundary site. Hence the system behaves like a 
time-homogeneous system where at each time a particle 
can exit with the effective rate P''^^ — 1/2, which corre- 
sponds to p'jI^ = 1 - f3^ff = 1/2 in TASEP. So we have 

lim^^oo Pi/'^ (^^) = 1/2 which we expect to be a good ap- 
proximation for all frequencies v ^ 1 much larger than 
the jump attempt rate. 

On the other hand, in the static case of zero fre- 
quency the system relaxes into the high density phase 
by a back-moving shock if the traffic light cycle starts 
with /3 — 0. This leaves the system with a bulk den- 
sity p = 1. If started with a green traffic light, (3=1, 
the system reaches the maximal-current phase with bulk 
density p = 1/2 [2l| . Continuity in frequency then gives 
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fyP (li) > 1/2 for all frequencies v provided that A > 1/2 
[27| . If the system is initially in the low density phase, 
then in the high frequency limit one has again an effec- 
tive right boundary density of 1/2, leaving the system 
in the low density phase with a bulk density A. In the 
zero- frequency case the limiting behavior depends again 
on how the period starts. If /3 = (red traffic light), the 
system fills up as described above and p = 1. On the 
other hand, for /? = 1, the system remains in the low 
density phase with a bulk density A. 

Notice that neither of these zero-frequency behaviors 
represents the zero-frequency limit i/ — s- shown in FigjSl 
For this limit we demonstrate in Sec. IIVI that 



andp^/-^H> p*. (2) 



This result is based on hydrodynamic limit arguments 
which are the subject of the following subsections. 



B. Hydrodynamic limit and mean field description 

The most basic question to be asked about the dy- 
namics of an interacting particle system is its large- 
scale behavior, i.e., how macroscopic equations of motion 
arise from its microscopic dynamics. By suitable coarse- 
graining of space and time the law of large numbers usu- 
ally guarantees that stochastic variables, in the present 
case the particle number in some interval (which under 
scaling becomes infinite on microscopic scale, but still in- 
finitesimal on macroscopic scale), turn into mean values 
whose temporal evolution satisfy some deterministic evo- 
lution equation in rescaled macroscopic time. Moreover, 
on macroscopic time scales the system is locally station- 
ary, i.e., all fast variables not captured in the evolution 
equation are locally stationary. This fact determines the 
precise form of the macroscopic equation, provided the 
stationary distributions are known. 

In the case of a conserved quantity the evolution equa- 
tion is a conservative pde of the form 



dp , d{j{p)) 



dr 



dx 



= 



(3) 



where a; is the rescaled space variable and r is rescaled 
macroscopic time [l^ [22] . The quantity j (p) is the 
particle current which on the macroscopic time scales 
for which ([3]) is valid takes its locally stationary value 
at density p{x,t). In the case of the ASEP one has 
= (p ~ 9)p(l ~ p)- Generally, for driven particle 
systems with a finite macroscopic current j{p) the ap- 
propriate hydrodynamic scale is the Euler scale t = ta 
X = ka, where a — > is the lattice constant. For the 
ASEP investigated here we get the inviscid Burgers equa- 
tion 



dp 

dr 



+ ip-q) 



djpjl-p)) 
dx 



(4) 



on the semiline a; < with traffic light boundary condi- 
tions on the right boundary x = 0. For the TASEP we 
set p = I and g = 0. 

Notice that in the present case of periodic boundary 
driving also the period has to be rescaled. The bound- 
ary conditions become p{Q,t) = i (1 + sign[sm(2Trflt)]) 
where = v/a is the rescaled frequency. I.e. the 
boundary stays open (p(0,rgreen) = ) during green 
light half-periods Tgreen C [0, |^], [T, ... and closed 
p(0, Tred) = 1 during the "red light" half-periods r^ed C 
[|^,T],[^,2T],.... Here we denote by T the complete 
period T 

The inviscid Burgers equation is the zero-viscosity 
limit of the viscous Burgers equation 



dp 
d^ 



dip{l-p)) ^d^P 



dx 



= D 



dx"^ 



(5) 



which can be solved in explicit form by a Hopf-Cole trans- 
formation for fairly general boundary conditions [1^ . For 
traffic light boundary conditions, however, such a solu- 
tion is difficult to obtain and we solve the coarse-grained 
time evolution of the ASEP by numerical integration. 
To this end we note that the exact microscopic operator 
equations of motion for the expected particle number Uk 
on site k read 

d 

Oli'^k) ^p{nk-i (1 - Uk)) - q{nk (1 - Uk-i)) 
- p{nk (1 - Uk+i)) + q{nk+i (1 - rik)) 

In this equation a one-point function (the expected den- 
sity {uk)) is coupled to two-point functions (on the r.h.s. 
of the equation), i.e. the equation is not closed. Writ- 
ing down an exact equation for the two-point functions 
introduces three-point functions and so on. This infinite 
hierarchy of equations is not directly tractable and some 
closure scheme must be employed for further analysis. 

In the mean field approximation for the ASEP, we 
neglect the correlations and approximate {rikUk+i) ~ 
{nk){nk+i) = PkPk+i etc. where pk = (rifc) is an av- 
erage particle density at site k. Using this approach for 
the exact microscopic evolution equation we obtain after 
some algebra 



d_ 
dt 



Pk = -{p- q) 



(1 - 2pk) 



Pk+l — Pk-1 



+ {p + q) 



Pk-i + Pk+i - '2pk 



(6) 



These equations are complemented with the traffic-light 
boundary conditions po {t) — pn (t) in ^ . 

Some comments are in order. Firstly, in the continuum 
limit one substitutes ka — x, with a ^ 1 being the lattice 
constant, e.g. pk+i (t) p{x + a,t). Taylor-expanding 
©, and using (p + q — 1) we obtain, after rescahng time 
as T = to, the Burgers equation 



dp 
d^ 



+ {p~q) 



d{p{l-p)) 
dx 



a d'^p 
2fe2 



(7) 
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Figure 6: Results of numerical integration of mean field equa- 
tions © for system size iV = 3000, = T"^ = 0.002, with 
initial condition p{x) = 0.5 after 10 full periods (thick curve). 
Thin curves show the subsequent density profiles p{x, t) after 
1, 2, 3, 4 and 5 full periods. The broken line shows the curve 
r PfT2)) exact in the limit ly ^ 0. 



Therefore, the mean field equation ([6]) can be viewed 
as a discretization of the viscous Burgers equation (O, 
with a constant discretization step a, and D = In 
fact, for numerical integration of the Burgers equation 
([5]) with the diffusion coefficient D, we shall be using the 
discretization scheme ^ with p = 1, q = 0, discretiza- 
tion step a — I and the coefficient 2D, instead of {p + q) 
in front of discrete second derivative. 

Secondly, for the case of weak hopping asymmetry 
lima^o i{p — q) = 1 one can obtain Eq. ([5]) from ^ by 
diffusive rescaling r = to? /D, x' = xa^/D in the hydro- 
dynamic limit a 0. In this case, the density does not 
evolve into shocks. A stationary travelling wave solution 
is a hyperbolic tangent with a step width proportional 
to the square root of the viscosity. Such smoothening 
of a shock may also be expected from the integration of 
the discrete mean field equation [B) In Fig. [5] we show 
the results obtained from numerical integration of the 
mean field equations ([6]) when a traffic light boundary is 
present at the origin. We see that the density profiles 
p{x, t) display a sawtooth structure which resembles the 
one observed for ASEP (see Fig. [3]). 

As a warning to readers not familiar with hydrody- 
namic scaling, we remark that the agreement between the 
mean field equation for the ASEP in the continuum limit 
and the rigorously derived Burgers equation is purely co- 
incidental. It arises from the fact that for the ASEP the 
stationary distribution has no correlations and hence the 
mean field equations become exact. In general lattice 
gases, e.g. in the KLS-model 0, one has correlations 
and a simple-minded mean field approach neglecting all 
correlations would produce a macroscopic equation that 
is in general not even qualitatively correct. 
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Figure 7: The universal curve F (solid curve) at the beginning 
of the green light period r = 0, T, 27, ...(thick curve) and 
for intermediate r values r/T = 1/4, 1/2,3/4 (broken, thick 
broken, thin curves, respectively). Inset shows F snapshots 
over larger scale. Note that microscopically the snapshots of 
average density profiles during the green light period pr — 
do not have discontinuity at the first site, but a boundary 
layer involving many sites, interpolating between p = 0.5 and 
PR — 0. However, the F snapshots show the limiting density 
profiles, rescaled by the period T, in the limit when T ^ oo 
[v = -fO). The boundary layer vanishes in this limit leading 
to a discontinuity at a; = 0. 



C. Stationary sawtooth structure 

The time-periodic sawtooth solution shown in FiglS]is a 
shock analog for the case of periodically changing bound- 
ary conditions which we now describe in detail. To this 
regard, we take the analytically tractable Burgers case 
(|3]) as a concrete example for explicit computation. The 
results for stationary (periodically repeating ) solutions 
obtained in this case are expected to be valid also for 
generic conservation law equation with convex j (p) . 

Looking at snapshots of density profiles p{x, t + mT) 
at times differing by multiples mT of a period, one 
notes that they fill some universal curve F, which has 
a characteristic sawtooth shape, see FiglS] The curve 
F-y(x, r) = r^(a;,T + T) "breathes" and returns to its 
original form after a full time period T, see Figl?] The 
index 7 denotes a phase < 7 < T at which the snap- 
shot of F is taken, with respect to the beginning of a 
green light interval. In the following we shall set 7 = 
and omit 7 for brevity of notation. 

Firstly, we describe the curve F and then prove its pe- 
riodicity in time. The curve consists of infinite number of 
sawteeth with heights decreasing away from the bound- 
ary. We denote the height, the base and the coordinate 
of /c-th sawtooth at the beginning of green light periods 
by and a;°, respectively, and the sawtooth angle 

by as shown in Fig[51 Each sawtooth is bounded by a 
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Figure 8: The universal curve F (solid curve) at the beginning 
of the green light period and variables gk,hk, otk used in the 
analysis (|9I12|I . 



shock discontinuity on the left and by a rarefaction wave 
on the right except the sawtooth fc = bounded on the 
right by a jam (caused by just finished red light period). 
Shock discontinuities move with the velocities given by 
Rankine-Hugoniot condition 



'^shock 

(k) 



i jhk) - j {9k) 

hk - gk 



(8) 



where j{p) is a flux function from ([3]). In the framework 
of stochastic driven systems j{p) is called the current- 
density relation or fundamental diagram. Shapes of rar- 
efaction waves are also determined by j {p). In the fol- 
lowing we shall consider a specific example ([4]). However 
we expect that our main results ()30|31|32p are applicable 
for arbitrary convex function j [p). 

For the specific case of inviscid Burgers equation, 



h^,g1, and for fc > are given by 




and for k = 



xl^kT[l-2K)\ll + - 



tan (al) = 2kT, 



7o" = ^ + Ac 



(9) 
(10) 

(11) 
(12) 

(13) 



where Ac = 
1/8. For k 

1-A, 



V2) /4 « 0.146 satisfies (1 - AJ = 
both hi 1- Xc + and w 

l^|r^ approach the limiting value 1 — Ac which is 



(2- 
> 1 



therefore the limit oi p {x,t) as x ^ —oo. This value may 
be identified with the effective time-independent bound- 
ary density p^/ . For rigorous definition of the boundary 
density see [131 • Here we only stress that the underly- 
ing particle system (ASEP) behaves as if it was joined 
at the origin with the reservoir of particles with the den- 
sity p'jj/ . Correspondingly, the average flux through the 
boundary is 



(jout) = j (pr) j (1 - Ac) 



(14) 



In the following we prove that T is a periodic function 
of time T{t) =T{T + t). In this respect it is sufficient to 
show that the heights and positions of all sawteeth will be 
the same after time T. However, since all sawteeth shocks 
discontinuities move to the left 9x„ (r) /dr < (this fol- 
lows from hn (t) > g„ (t) > 1/2 and ([5])), sawteeth can- 
not return to their original places. Instead, after time T a 
sawtooth k will take the place of a former sawtooth k+1, 



I.e. 



hk (T) 



'''k+l 



9k (T) = 5^+1 and Xk (T) 



for all k, see FigUl At the boundary, the structure with 
a jam at position Xq must be regenerated after complete 
green hght r C [0, T/2] and red right r C [r/2, T] period. 
To proceed, note that the inviscid Burgers equation (|4]) 
has two basic solution types: (a) between two consecu- 
tive homogeneous states p_ < p+ a shock discontinuity is 
formed moving with velocity v shock — P-^ P+ and (b) 
between two consecutive homogeneous states p_ > a 
rarefaction wave, p{x,t) = p{x/t) — 1/2 — x/{2t) is 
formed. This information is enough to predict the evolu- 
tion of sawteeth structure in Fig[51 consisting of shocks 
and rarefaction waves. In particular, the shock velocity 
of a k-th sawtooth 



dxk 
'd7 



and 



dr 



1 - hk (t) - gk (t) 



tan(Q;fe (r)) 



(15) 



(16) 



In the following we shall explicitly indicate time- 
dependent quantities e.g. hk (r) , ak (t) while h^, will 
denote their initial values hi = hk (0) etc., given by 
(dnil-dni). The time interval < r < T will be con- 
sidered. For fc > Eqs. (|16ll2p are trivially solved: 
tan ak (t) — 2 (t + kT). For hk (t) , gk (t) we can write 
the equations 



-Xk (t) 



hk (t) - 



Y = tan ak (r) = 2 (t + fcT) (17) 



-Xk (t) 

1 

2 



gk (t) 



tanafc+i (r) = 2 (r + (fc + 1) T) . 
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Multiplying pT)) by hk (t) — i, differentiating with re- 
spect to T and using (|15p we obtain 

2 (r + kT) + /ifc (r) - (r) = 0. 

Analogously, for gk (t) we get an ordinary differential 
equation 

2{t + {k + 1)T)^^^ + gk{T) - hk (r) =0 

OT 

with the initial conditions gk{0) — g^, hk{0) = h^- 
These equations can be integrated to give 



(18) 



1 



)T + T ^ 

In particular, at time t = T 

hk {T) = hi + {hi - gl) yk{k + 2)-k- l} (19) 



gk{T)=gl + {hl~gl)k 



k + l 



y^k {k + 2) 



- 1 



(20) 



It can be verified that (T) = ctl+i^ ^fe (-^) ^ 
and gk (T) — g^^i, thus ensuring the regeneration of the 
curve r after a period T. Special attention should be 
given to the boundary region A; = 0. The equation for 
50 (■'■) reads 



'Xo (r) 



90 (r) 



- = 2(T + t) 



Multiplying by denominator and differentiating with 
respect to t we obtain 2 (t + T) ^^§^ + 2 (go (r) - i) = 
_dxo^^ At point xq (t) there is a jump, 
P (xo (t) - 0, t) = go (r) and p(xo (r) + 0, r) = 1 
at the right, consequently 

^ = -^"^^l^^-.Y"" - -.90 (r). Substituting, we ob- 
tain 

2(r + r)^+,o(r)-l = 0. 

OT 

Solving the latter with the initial condition gQ (0) = g^, 
we get 



50 (r) = 1 + (gS - 1) 



T + r' 



(21) 



describing the shock propagation during < r < T. 
Alongside, the jam joining the boundary will start to 



dissolve by rarefaction wave p{x, r) 



^ during the 



green light period < r < -j . During the red light period 
Y < T < T new jam appears at the boundary a; = 
and propagates inside. We shall denote its coordinate by 
XQ (r), xg (^) = 0. The base of the new jam, denoted 
by G(T),G(|)-i will obey 

-XG (r) 



G(r) 



= 2r. 



Using dxG{T)/dT = ^G{t), we obtain an equation 
2t^ + G (r) - 1 = 0, solved by G(r) = l-i^^ 
for f < T < T. At time r = T one has G (T) = 1- ^ = 
i + Ac = 5o, restoring the initial shape at r = 0, see (|13p . 
g^ is determined from the solution (|2ip . go (2^) = .9i- At 
time T = T the jam, which was initially at position Xq, 
reaches point x^, consequently h^ = 1. The remaining 
relations ([9l)- pT|) are obtained recurrently from (|19l20p . 
Thus the proof of periodicity of the sawtooth structure 
is completed. 



D. Steady state selection 

In the previous subsection we have proved stationarity 
of the saw tooth state under periodic driving, but we did 
not address the question whether this stationary state 
is actually reached for any initial state characterized by 
the initial density A. In order to investigate this problem 
of steady state selection we consider an initial state con- 
sisting of homogeneous state on the left p_ {x) = A and 
sawtooth structure on the right, joined by a shock, and 
demonstrate that it is analogous to a shock between the 
two homogeneous states P-{x) = A and p+ (a;) = 1 — Ac. 

The latter shock moves with the velocity v — 
(j (A) - j (1 - Ac)) / (1 - A - Ac) = Ac-A, and wiU travel 
a distance (Ac — A) T after time T. In particular, the 
shock will be stationary for A = Ac. Let us prove this 
feature for the shock between the homogeneous state 
P-{x) = Ac, for X < Xn and sawtooth structure with 
n complete sawteeth p+{x) — F, for Xn < x < 0. Heights 
of all sawteeth hk, k < n will satisfy pT|) . but the velocity 
for the n-th sawtooth satisfies 



dxk _ hk{T){l - hk{T)) - Ac(l - Ac) 



= l-hk (t) - Ac 



Ot hk{T) - \c 

(22) 

Multiplying ([TT]) for fc = n by /i„ (r) — i, differen- 
tiating with respect to r and using (I22p we obtain 
2 (t + nT) ^^gj^-* + hn (r) — A — 0, whose solution is 

hn (r) = Ac + - Ac) \J nT+T ' 'o^deA to prove that 
the shock will return to the original position one after 
time T, i.e. /J" %f-dT = 0. Substituting the solution for 
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hn (t) into ([22|l . and integrating over the period, we have 
r ^dr = (2Ae-l)r+2T«-Ae)(Vn(" + = 0. 

Jo OT 

(23) 

Substitution ([51) satisfies the above equation. 

Analogously, one proves that the shock between A 7^ Ac 
and the sawtooth structure, after period T, will advance 
if A > Ac or retreat towards the boundary if A < Ac- 
The shock position after time T is determined by a bal- 
ance equation. Namely, the extra mass gained by the 
shock, AA/ = (p (a;, T) — p (x, 0)) dx, is equal to 
the difference between the ingoing and outgoing currents 
(jm-Jo«t)T=(i(A)-j(l-Ac))r, 

AM = (j,„ - jout) T = (j (A) -3 {I- Ac)) T. 

The explicit calculations are not illuminating and are 
omitted for brevity. As a guide to the eye, see FiglH] 
The most rapid progression of the sawtooth structure is 
achieved for jin = jmax = 1/4 for which the number of 
sawteeth n increases roughly by two each five cycles, see 
FigH 

The variations of the boundary density considered in 
this paper so far (Eq.(IT])) followed a square wave in which 
the signal is green for a fraction / = Tgreen/T = ^ of the 
period. What happens if we vary the fraction / between 

and 1? While the cases / = and / = 1 are obvious and 
were discussed before Eq®, the general case < / < 1 
produces a sawteeth structure r(/), qualitatively similar 
to the one for f — 1/2, compare Figs l8l6l and FigO The 
limiting curve r(/) converges in the bulk to the value 

1 — A (/) = (1 + Vl — /) /2, determined by the outgoing 
flux in the vanishing frequency limit (jo«t(/)) = /jmax, 
see the discussion after Eg ([50)1 . Note that A (/) satisfies 
A (/) (1 - A (/)) = {joutif)) = //4. Analytic analysis for 
arbitrary / can be carried out analogously to the / = 1/2 
case. In particular, the heights of the sawteeth (/) 
are determined by the solution of the Eq. (P5)) where Ac 
is substituted by A (/) = (l - V^~~f) /2, 

+ (24) 

valid for n > such that ft," (/) < 1. This is always the 
case if / > 1/2 ( duration of green signal is larger then 
the duration of the red signal). In the opposite case, 
/ < 1/2, additional plateau appear with the saturated 
density p = 1 on the curve r(/), apart from the very 
first one caused by the red traffic light at the boundary, 
see upper curve in FiglH) A thorough analysis of these 
more complicated structures is beyond the scope of the 
present paper. As a guide for an eye, see FiglQl where 
the curves T (/) for / = 0.25, 0.5, 0.65 are shown. 

Let us stress once more the universality aspect of the 
limiting curve F shown in Figs lTIBI The shape of the 
curve is independent on T provided that T is sufficiently 
large = 1/ — > 0. It also does not depend also on 
ASEP bulk rates p, q provided that p > q (drive towards 




Figure 9: Sawteeth density profiles for different fractions of 
green signal / = Tgreen/T = 0.35, 0.65 Thick (/ — 0.35) and 
thin (/ = 0.65) curve show respective density profiles p{x,t) 
from mean field equations (|6} evolving from the initial con- 
dition p{x, 0) — 0.5 after 10 full periods. In MP calculations, 
1/ = T^^ = 0.001. Broken lines show the limiting curves F (/), 
exact in the limit u 0. The sawteeth heights are given by 
(|24p . All the curves are shown at the moment of time when 
the red light turns off. 

the right boundary), since F is described by the equation 
(HI). Qualitatively, F does not depend on relative dura- 
tion of the green light period with respect to the whole 
period T (equal to 1/2 in the present study), as argued in 
the previous paragraph. Moreover, for other models with 
the convex current-density relation j{p) and traffic light 
boundary conditions we expect the existence of a simi- 
lar curve with sawtooth structure, with model-dependent 
shape of the sawteeth. The distance between nearest saw- 
teeth is determined by j{p), see (I5T|) . The sawteeth curve 
in the bulk will converge to a value p^'^ (v), determined 
by averaged outgoing flux through (jout) = j [p]^/ [v]^ ■ 
The outgoing flux in the limit — > will be given by 

dsni). 



IV. FINITE SIZE ASEP WITH TRAFFIC LIGHT 
BOUNDARY CONDITIONS AT ONE OR BOTH 
BOUNDARIES. 

Here we discuss stationary behaviour of sufficiently 
large but finite system of size N with open boundaries. 
It is intuitively clear that the periodically changing con- 
ditions at a boundary will generate a sawteeth structure 

with a typical sawtooth size A(i/) w 1- 

We shall consider the case 3> A. If A^ < A, then one 
boundary will influence another one during a periodic 
cycle. 

Firstly, recall the well-known results for the TASEP 
model with time-independent boundary rates, solved ex- 
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Figure 10: Stationary density as function of the frequency v 
for fixed left boundary density pL = 0.3, from Monte Carlo 
simulations (diamonds with errorbars), and from integration 
of Burgers equation (0 with D = 2, 0.9 (solid lines) and D — 
0.5 (broken line). The dotted line shows (patat) computed 
from stationary flux measurements. Fluctuations at small 
frequencies are due to flnite size effects. 



actly in [13|,llJ|. In the TASEP, a particle can be injected 
at first site k = 1 from the left boundary reservoir with 
the rate a and be extracted from the last site k = N 
with the rate /3. This corresponds to coupling at the left 
with the reservoir of particles with density pL = a and 
on the right with the reservoir of particles with the den- 
sity PR = 1 — (3. In the range of parameters < a, /? < 1 
the stationary states are characterized by average ho- 
mogeneous particle distribution in the bulk with density 
Pstat (pl, Pr)- The stationary densities obey an extremal 
principle for the stationary flux |5j] , 

„ f min[p^^p^] j (p) if PL < PB. 

max[,,^,^]j(p) \ipL>PR ' ^^^^ 

which, being applied for the case of ASEP j {p) — 
p{l — p), yields three different phases. 

Low density (LD) pstat = PL, for pL ^ a < 1/2, pB. < Pl 

(26) 

High density (LD) p^tat = PR, for pR > ^f^^PR > Pl 

(27) 

Max. current (MC) pstat = 1/2, for pL > 1/2, pR < 1/2. 

(28) 

In the case of fixed left boundary conditions pL — const 
and periodically changing right boundary conditions (e.g. 
traffic light boundary conditions) the pR in (|25ll28p has 
to be replaced by the effective boundary density p^-^ (z^) 
and j Stat, Pstat by a time-averaged flux and density in 
a stationary state {jstat),{Pstat} respectively. Note that 



Figure 11: Phase diagram of ASEP with open boundaries, il- 
lustrating (|26I28|) . Solid line indicates discontinuous LD/HD 
transition. Dashed line borders the region reachable by ap- 
plying the traffic light boundary conditions with arbitrary 
frequencies at both boundaries. 



since p^J^ {v) > 1/2 for any i', only LD and HD phases 
can be observed. 

In particular, (|25ll28p with the latter substitution pre- 
dicts a discontinuous change of stationary density from 
{Pstat) =PL< 1/2 (LD phase) to (pstat) = Pr" M (HD 
phase) at a transition point 1 — pL — p'r^ {v)- Indeed, 
keeping p^ fixed and changing v, one observes this phase 
transition at the predicted point, see FigfTOl 

Analogously, keeping fixed right boundary conditions 
PR = const and applying traffic light boundary condi- 
tions at the left boundary, one has to define the effec- 
tive left boundary density p'j^^ (v). The latter, due to 
particle-hole symmetry of the TASEP, is given simply by 



(29) 



varies between 



Consecutively, the 

Pl^ +0) = to p^" (i/ ^ oo) = 1/2. 

Finally, applying traffic light boundary conditions at 
both boundaries with the frequencies v and v' on the left 
and on the right, respectively, effective boundary reser- 
voirs p']^^ {v) at the left and p''^^ {v') at the right are 
created. Again, one finds the phase diagram applying 
the rule (P5|) with the replacements pL p^l^ {v),Pr 
Pr/ (^')- Taking into account ([29]) . one predicts LD ((26|) 
phase (pstat) = p'^R^ {v) for v > v' and HD (|^ phase 
{Pstat) = p'^L^ {v) = I - p^l^ {v) for V < v'. Due to 
the range of variance of the effective boundary densi- 
ties p*^/^ {v) < 1/2, 1/2 < p^-^ (i/), the maximal current 
phase ([28|) cannot be reached except at one point, see 
FigHH 

In addition to Monte-Carlo simulations, we integrated 
numerically mean field equations In the limit of in- 
finitely small frequencies — > the solution of the mean 



10 



field equations converges to the solution of the invis- 
cid Burgers equation, see Figs l6l9l For finite frequencies 
J/ > 0, we observe qualitative agreement between mean 
field and Monte-Carlo density profiles, see FigHl How- 
ever, the mean field solution apparently fails to to predict 
the exact location of the phase transition frequency, see 
the curve with D — 0.5, marked MF, in FigdOl Indeed, 
as argued after Eq. ([7|) , meanfield equation can be viewed 
as a discretization of viscous Burgers equation with the 
diffusion coefficient D = 0.5. The mean field approach 
fails quantitatively because it neglects correlations which 
are present in the sawtooth state. On the other hand, 
numerical integration of the Burgers equation ([5]), keep- 
ing the diffusion coefficient D a free parameter, shows 
that the the effective boundary density (and consequently 
the critical frequency) depends on D. Manipulating Z?, 
one can obtain a better agreement with the Monte Carlo 
simulations, see FigfTOl It might seem from the FigfTOl 
that the choice D = 0.9 fits the Monte-Carlo data well. 
However, the deviation between Monte-Carlo data and 
Burgers equation with D = 0.9 for small frequencies ly 
is substantial and it can be seen by comparing graphs of 
type Figlinifor different (not shown). 



V. CONCLUSIONS 

We have provided a hydrodynamic description of the 
semi-infinite ASEP with traffic light boundary condi- 
tions. We find a time-periodic stationary sawtooth solu- 
tion which is described in detail. We have also addressed 
the question of steady state selection, starting from some 
initial density A. The picture that emerges is similar to 
that of the usual ASEP with constant effective reservoir 
density that we have determined. Despite the sawtooth 
structure of the solution, the time averaged density is at 
sufficiently large distance from the boundary is given by 
the p{x,t) = p'^-^-^ , both in low and high density regime. 
Indeed, considering the motion of the shock as an effec- 
tive the single-particle problem in an external potential 
[13, UHl this observation is reminiscent of the motion of a 
Brownian particle in a periodically driven stochastic sys- 
tem [26|. Our result shows that effective potentials may 
arise also in interacting many-body systems as a result 
of periodic driving. 

Our derivation is based on the ASEP as a specific 
example, but remains valid for generic driven diffusive 
systems with convex current-density relation. There 
are several quantitative conclusions that one can draw 
from the exact hydrodynamic treatment presented above. 
Firstly, note that the average outgoing flux, see ((M]) 
(jout) = 1/8 is two times smaller than the maximal flux 
maxj (p) = jinax = 1/4. The relation 



(jout) = 



Jn 



(30) 



is not a casual, but a rather generic one: during the 
red light periods, {joutjred = 0, and an extended jam 



at the boundary forms. During the green light period, 
{jout) green = jmax, bccausc the outflow from a jam 
is governed by a maximization principle {jout) green = 



[04] j iP) = jn 



Per full period, one obtains 



Secondly, the distance between neighbouring sawteeth 
rapidly approaches a constant, \x}^+i — Xk \ ~ (1 — 2Ac) T 
for fc 3> 1, as follows from (fTT|) . The value of the con- 
stant has simple physical origin: maximum and min- 
imums of sawtooth structure F, amd gu approach 
the effective boundary density value hk,gk ~ 1 — Ac ± 
0{l/k) — Pa ± 0{l/k). Hence, the velocities of the 
discontinuities for large k approach the group veloc- 

Ity Vgroup (PR) - limfc^oo ht{T}-gk{T} ^ 



Xk+l ^ Xk - 



By periodicity r(0) = F (T) requirement 

'group (Pr) T, or 



\Xk+l - Xk\ 

T 



(31) 



The distance between sawteeth converges monotoni- 
cally and rapidly to the predicted value: indeed, as 
follows from (fTTj) . liuik^ao \xk — xo\/T = {k + ^ — 



2Ac 



:) (dj/dp)^^^ 



Consequently, the sum of all de- 



viations does not exceed 10% of the predicted limiting 
distance A. Hence, one can measure derivative of the 
flux {dj/dp)p^p^ directly by measuring the distance be- 
tween the sawteeth. Estimate of {dj/dp)^^^^ from ([3T|) 
for A; = 1 (the flrst and the best-visible sawtooth) induces 
relative error less than 4%. 

Finally, the amplitude of density variations 5p (x) can 
be estimated as — 3^1 at a distance \xk\ ~ TfcA. From 
®,((Tni) we obtain 



5pi-TkA) 



l-2Ac 
2k 



O 



(32) 



In most of the paper, a semi-infinite system was consid- 
ered. The analogy of the sawtooth structure with a shock 
and the picture of steady state selection that emerges 
allows us to consider also finite systems with two open 
boundaries through which particles can enter or leave the 
system. We argue that the traffic light boundary condi- 
tion represents a domain of the full phase diagram that 
includes part of the first-order transition between low- 
and high-density phase. The maximal-current phase is 
reached only in the point where it meets the end of the 
first order transition line. 
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